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METHODS FOR TREATING BIPOLAR MOOD DISORDER 
ASSOCIATED WITH MARKERS ON CHROMOSOME 18q 
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INTRODUCTION 

Background 

Bipolar Mood Disorder (BP) 

Manic-depressive illness, or bipolar mood disorder (BP), is characterized 
by episodes of elevated mood (mania) and depression and is among the most 
prevalent and potentially devastating of psychiatric syndromes. The most severe 
and clinically distinctive forms of BP are BP-I (severe bipolar mood disorder) and 
SAD-M (schizoaffective disorder manic type), and are characterized by at least 
one full episode of mania, with or without episodes of major depression (defined 
by lowered mood, or depression, with associated disturbances in rhythmic 
behaviors such as sleeping, eating, and sexual activity). A milder form of BP is 
BP-IL bipolar mood disorder with hypomania and major depression. BP-1 often 
co-segregates in families with more etiologically heterogeneous syndromes, such as 
unipolar major depressive disorder (MDD), which is a more broadly defined 
phenotype. See Mclnnes, L.A. and Freimer, N.B., Mapping genes for psychiatric 
disorders and behavioral traits, Curr. Opin. in Genet, and Develop., 5:376-381 
(1995). 
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Treatment of Individuals With Bipolar Mood Disorder 
An estimated 2-3 million people in the United States are affected by BP-L 
Currently, individuals are typically evaluated for bipolar mood disorder using the 
clinical criteria set forth in the most current version of the American Psychiatric 
Association's Diagnostic and Statistical Manual of Mental Disorders (DSM). 
Many drugs have been used to treat individuals diagnosed with bipolar mood 
disorder, including lithium salts, carbamazepine and valproic acid. However, 
none of the currently available drugs is able to treat every individual diagnosed 
with severe BP-I (termed BP-I) and drug treatments are effective in only 
approximately 60-70% of individuals diagnosed with BP-I. Moreover, it is 
currently impossible to predict which drug treatments will be effective in particular 
BP-I affected individuals. Commonly, upon diagnosis affected individuals are 
prescribed one drug after another until one is found to be effective. Early 
prescription of an effective drug treatment is critical for several reasons, including 
the avoidance of extremely dangerous manic episodes and the risk of progressive 
deterioration if effective treatments are not found. Also, appropriate treatment 
may prevent depressive episodes in BP-I individuals; these episodes are also 
dangerous and are characterized by a high suicide rate. The high prevalence of the 
disorder, together with frequent occurrence of hospitalizations, psychosocial 
impairment, suicide and substance abuse, has made BP-I a major public health 
concern. 

Genetic Basis for Bipolar Mood Disorder 

Mapping genes for common diseases believed to be caused by multiple 
genes, such as BP-I, may be complicated by the typically imprecise definition of 
phenotypes, by etiologic heterogeneity, and by uncertainty about the mode of 
genetic transmission of the disease trait. With psychiatric disorders there is even 
greater ambiguity in distinguishing between individuals who likely carry an 
affected genotype from those who are genetically unaffected. For example, one 
can define an affected phenotype for BP by including one or more of the broad 
grouping of diagnostic classifications that constitute the mood disorders: BP-I, 
SAD-M, MDD, and BP-II. 
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Thus, one of the greatest difficulties facing psychiatric geneticists is 
uncertainty regarding the validity of phenotype designations, since clinical 
diagnoses are based solely on clinical observation and subjective reports. Also, 
with complex traits such as psychiatric disorders, it is difficult to genetically map 
the trait-causing genes because: (1) the BP-I phenotype doesn't exhibit classic 
Mendel ian recessive or dominant inheritance patterns attributable to a single 
genetic locus, (2) there may be incomplete penetrance i.e., individuals who inherit 
a predisposing allele may not manifest the disease; (3) the phenocopy phenomenon 
may occur, i.e., individuals who do not inherit a predisposing allele may 
nevertheless develop the disease due to environmental or random causes; (4) 
genetic heterogeneity may exist, in which case mutations in any one of several 
genes may result in identical phenotypes. 

The existence of one or more major genes associated with BP-I and with a 
clinically similar diagnostic category, SAD-M (schizoaffective disorder manic 
subtype), is supported by segregation analyses and twin studies (Bertelson et al., 
1977; Freimer and Reus, 1992; Pauls et al.. 1992). However, efforts to identify 
the chromosomal location of BP-I genes have yielded disappointing results in that 
reports of linkage between BP-I and markers on chrci..osomes X and 1 1 could not 
be independently replicated nor confirmed in the re-analyses of the original 
pedigrees (Baron et al., 1987; Egeland et al., 1987; Kelsoe et al., 1989; Baron 
et al.. 1993). Recent investigations have suggested possible localization of BP 
genes on chromosomes 18 (pericentromeric region) and 21 q, but in both cases the 
proposed candidate region is not well defined and there is equivocal support for 
either location (Berrettini et al. (1994) Proc. Natl. Acad. Sci. USA, 91, 5918- 
5921, Murray. J.C., et al. (1994) Science 265, 2049-2054; Pauls et al.. Am. J. 
Hum. Genet. 57:636-643 (1995): Maier et al.. Psych. Res. 59:7-15 (1995); Straub 
et al., Nature Genet., 8:291-296 (1994)). 

Despite abundant evidence that BP has a major genetic component, linkage 
studies have not yet succeeded in definitively localizing a BP gene. This is mainly 
because mapping studies of psychiatric disorders have generally been conducted 
under a paradigm appropriate for mapping genes for simple Mendelian disorders, 
namely, using linkage analysis in the expectation of finding high lod scores that 
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definitively signpost the location of disease genes. The follow up to early BP 
linkage studies, however, showed that even extremely high lod scores at a single 
location can be false positives. See Egeland, et al., Nature 325:783-787 (1987); 
Baron et al.. Nature 326:289-292 (1987); Kelsoe et al., Nature, 342:238-243 
(1989); and Baron et al.. Nature Genet. 3:49-55 (1993). These earlier studies 
used largely uninformati\ c markers and did not use stringent criteria for 
identifying affected individuals. 

Linkage Disequilibrium Analysis 

Linkage disequilibrium (LD) analysis is a powerful tool for mapping 
disease genes and may be particularly useful for investigating complex traits. LD 
mapping is based on the following expectations: for any two members of a 
population, it is expected that recombination events occurring over several 
generations will have shuffled their genomes, so that they share little in common 
with their ancestors. However, if these individuals are affected with a disease 
inherited from a common ancestor, the gene responsible for the disease and the 
markers that immediately surround it will likely be inherited without change, or 
IBD ("identical by descent"), from that ancestor. The size of the regions that 
remain shared (i.e. IBD) are inversely proportional to the number of generations 
separating the affected individuals and their common ancestor. Thus, "old" 
populations are suitable for fine scale mapping and recently founded ones are 
appropriate for using LD to roughly localize disease genes. (Houwen et al., 1994, 
in particular FIG. 3 and accompanying text). Because isolated populations have 
typically had a small number of founders, they are particularly suitable for LD 
approaches, as indicated by several successful LD studies conducted in Finland (de 
la Chapelle, 1993). 

LD analysis has been used in several positional cloning efforts (Kerem et 
al., 1989; MacDonald et al., 1992; Petrukhin et al., 1993; Hastbacka et al., 1992 
and 1994), but in each case the initial localization had been achieved using 
conventional linkage methods. Positional cloning is the isolation of a gene solely 
on the basis of its chromosomal location, without regard to its biochemical 
function. Lander and Botstein (1986) proposed that LD mapping could be used to 
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screen the human genome for disease loci, without conventional linkage analyses. 
This approach was not practical until a set of mapped markers covering the 
genome became available (Weissenbach et al.. 1992). The feasibility of genome 
screening using LD mapping is now demonstrated by the applicants. 

Identification of the chromosomal location of a gene responsible for causing 
severe bipolar mood disorder can facilitate diagnosis, treatment and genetic 
counseling of individuals in affected families. 

Due to the severity of the disorder and the limitations of a purely 
phenotypic diagnosis of BP-I, there is a tremendous need to genetically subtype 
individuals with BP-I to confirm clirucal diagnoses and to determine appropriate 
therapies based on their genotypic subtype. 

The present invention comprises using genetic linkage and haplotype 
analysis to identify an individual having a bipolar mood disorder gene on the long 
arm of chromosome 18 (18q). In addition, the present mvention provides markers 
linked to a gene responsible for susceptibility to bipolar mood disorder that will 
enable researchers to focus future analysis on that small chromosomal reg.on and 
will accelerate the sequencing of a bipolar mood disorder gene located at 18q. 



20 SUMMARY O F THF 1NVFMTTOM 

The present invention is directed to methods of detecting the presence of a 
bipolar mood disorder susceptibility locus in an individual, comprising analyzing a 
sample of DNA for the presence of a DNA polymorphism on the long arm of 
chromosome 18 between the markers D18S469 and D18S554, wherein the DNA 
polymorphism is associated with a form of bipolar mood disorder. The invention 
includes the use of genetic markers in the roughly 6-7 cM region between the 
markers D18S469 and D18S554 on 18q. inclusive, to genetically diagnose bipolar 
mood disorder in individuals and to confirm phenotypic diagnoses of bipolar mood 
disorder. 
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In a further embodiment, the invention provides methods of classifying 
subtypes of bipolar mood disorder by identifying one of more DNA 
polymorphisms located within the region between markers D18S469 and D18S554. 
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inclusive, and analyzing DNA samples from individuals phenotypically diagnosed 
with bipolar mood disorder for the presence or absence of one or more of said 
DNA polymorphisms. 

In yet a further embodiment, the methods of the invention include a method 
5 of treating an individual diagnosed with bipolar mood disorder comprising 

identifying one or more DNA polymorphisms located within 18q region between 
markers D18S469 and D18S554, and analyzing DNA samples from individuals 
phenotypically diagnosed with bipolar mood disorder for the presence or absence 
of one or more of the DNA polymorphisms, and selecting a treatment plan that is 
10 most effective for individuals having a particular genotype within the 6-7 cM 

region of chromosome 18. 



BRIEF DESCRIPTION OF THE DRAWINGS 
15 FIG. 1 is a pedigree chart showing two families, CR001 and CR004. 

Affected individuals are denoted by black symbols, deceased individuals by a 
diagonal slash. A schematic of each individual's haplotype (where available) is 
shown below the ID number. Recombinations are denoted by n -x"; 
consanguineous marriages by a double bar, and the conserved haplotype as black 
20 shading within the haplotype bars. The larger conserved region for CR004 is 

stippled, the larger conserved region for CR001 is indi. ated by a dashed outline. 
A "I" underneath the haplotype bars indicates inferred haplotype. A "?" indicates 
phase is uncertain. The connection between CR001 and CR004, dating to an 18th 
Century founding couple, is indicated by the dashed lines joining individuals III-6 
25 and 1-4. 

FIG. 2 is a table of lod scores for markers covering the entire human 
genome that exceeded the arbitrary coverage thresholds. Lod scores are shown for 
markers on chromosome 18, including one marker in the vicinity of the 18q22-q23 
30 region: D18S1161. 
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FIG. 3 depicts the extent of marker coverage used in the pedigree genome 
screening study for each chromosome. Coverage is defined as regions for which a 
lod score of at least 1 .6 would have been detected (in the combined data set) for 
markers truly linked to BP-I under the model employed. Areas that remain 
uncovered (at this threshold) are unshaded. Markers for which lod scores were 
obtained that exceeded the empirically determined coverage thresholds in CR001, 
CR004. or the combined data set, are shown at their approximate chromosomal 
location. The symbols to the right of the chromosome indicate the thresholds 
exceeded at that marker: a circle signifies that the lod score at a marker exceeded 
the threshold of 0.8 in CR001, a diamond signifies that the lod score exceeded the 
threshold of 1.2 in CR004, and a star signifies that the lod score exceeded the 
threshold of 1.6 in the combined data set. 

FIGS. 4A and 4B depicts the Lod score for the maximum likelihood 
estimate of theta in the combined sample for the 473 microsatellite markers typed 
in the pedigree genome screen. The MLEs of theta are represented by different 
colors: red = theta < 0.10; green = 0.10 < theta < 0.40; blue = theta 
>0.40. Note that the scale for the x-axis (distance from pter) changes with 
chromosomes. 

FIG. 5A: Affected members of two families CR001 and CR004 with 
depiction of the shared marker haplotypes observed. The unshaded area outlined 
in solid lines is the more extensive haplotype conserved in CR004, the unshaded 
area outlined in dashed lines is the more extensive haplotype conserved in CR001. 
The shaded area indicates a haplotype that extends from D18S1009 to D18S554 
and is common to both larger haplotypes. The ID numbers in the first column 
refer to the pedigree displayed in FIG. 1. In the remaining columns are the allele 
sizes at the indicated markers. indicates an uncertain haplotype, indicates 
that individual received two copies of the shaded haplotype (both haplotypes are 
displayed); "@" indicates an inferred haplotype. The markers used for 
haplotyping are shown at the top of the figure with inter-marker distances in cM. 
The marker order towards qter, is: D18S64. D18S55, D18S61. D18S485, 
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D18S870, D18S469, D18S1161. D18S1121, D18S1009, D18S380, D18S554, 
D18S462, D18S461, D18S70. 

FIG. SB: The other haplotypes that could be unambiguously reconstructed 
in the pedigree in FIG. 1 are shown. Those that are present in affected individuals 
are indicated by the ID numbers at the left of the diagram. "NT" indicates that an 
individual was not typed for a given marker. 

FIGS. 6A and 6B depicts publicly available markers in the 18q22-q23 
region of chromosome 18 along with their locations on a genetic linkage map and 
the radiation hybrid map. 

DESCRIPTION OF SPECIFIC EMBODIMENTS 
The recent availability of highly polymorphic, genetically mapped markers 
covering the human genome (Weissenbach, J., et al. (1996) Nature 359, 794-801, 
Murray, J.C., et al. (1994) Science 265, 2049-2054, Gyapay, G., et al. (1994) 
Nature Genet 7,246-339) has enabled the development of a multi-stage paradigm 
for mapping genes for complex traits. In the first stages, complete genome 
screening (e.g. through lod score analysis) is used to identify possible localizations 
for disease genes. Subsequently, the regions highlighted by the screening study 
are more intensively investigated to confirm the initial localizations and delineate 
clear candidate regions. Finally, fine mapping methods (such as haplotype or 
linkage disequilibrium (LD) analysis) or candidate gene approaches are used for 
positional cloning of disease genes. 

Our genome screening study for BP employed the following strategies. 
Unlike previous genetic studies of BP, only those individuals with the most severe 
and clinically distinctive forms of BP (BP-I and schizoaffective disorder manic 
type, SAD-M were considered as affected, rather than including those diagnosed 
with a milder form of BP (BP-II) or with unipolar major depressive disorder 
(MDD). Two large pedigrees (CR001 and CR004) were selected from a 
genetically homogeneous population, that of the Central Valley of Costa Rica (as 
described in Escamilla, M.A., et al., (1996) Neuropsychiat. Genet. 67, 244-253, 
and in Freimer, N.B., et al. (1996) Neuropsychiat. Genet, 67, 254-263, both 
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incorporated by reference herein). The entire human genome was screened for 
linkage using mapped microsatellite markers and a model for genetic analysis in 
which most of the linkage information derived from affected individuals. The goal 
of this stringent linkage analysis was to identify all regions potentially harboring 
5 major genes for BP-I in the study population. Empirically determined lod score 

thresholds (using linka-e simulation analyses) were derived, to suggest regions 
worthy of further investigation. 

Identification of all suggestive regions and weighing the relative importance 
of findings required complete screening of the genome. The coverage approach 
10 was developed to gauge the progress of this effort. Conventionally, the 

thoroughness of genome screening is evaluated by excluding genome regions from 
linkage under given genetic models. This approach, which is highly sensitive to 
misspecification of genetic models, maybe poorly suited for genome screening 
studies of complex traits; it is tied to the expectation of finding linkage at a single 
15 locus and demonstrating absence of linkage at all other locations in the genome. 

Additionally, exclusion analyses do not differentiate between genome regions 
where linkage is not excluded because markers are uninformative in the study 
population from those in which the genotype data are simply ambiguous. In 
contrast, the coverage approach is designed for studies aimed at genome screening 
20 rather than for studies whose goal is to demonstrate a single unequivocal linkage 

finding, and it provides explicit data regarding the informativeness of markers in 
the study pedigrees. Its use lessens the possibility that one would prematurely 
dismiss a given genome region as being unpromising for further study. 

Because the exact genetic length of chromosomes is not clearly established, 
it is impossible to be certain that one has screened the entire genome. Although we 
report coverage of about 94% of the genome (under the 90%) dominant model) at 
the thresholds described above, this probably represents an underestimate. The 
remaining coverage gaps in our study occur predominantly at or near telomeres; as 
the upper bound estimates for the length of each chromosome were used, it is 
likely that the actual coverage gaps in these regions are smaller than or 
conservative assessment. 
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The presence of consistently positive lod scores over a given region was 
considered to be of greater significance than isolated peak lod scores. Such 
clustering suggests true co-segregation of markers and phenotypes (i.e. alleles are 
shared identical by descent rather than identical by state) and is more readily 
observed in analyses of a few large pedigrees (as in our study) than in examination 
of several smaller families. The data presented herein indicates clustering of 
positive lod scores in the region of the 18q22-q23. The screening data provide 
obvious indication of a BP-I susceptibility locus in 18q22-q23 as suggestive lod 
(logarithm of odds) scores were obtained there over a region of about 40 cM. 

The genome screen was conducted in two stages. The Stage I screen 
identified areas suggestive of linkage, so that those areas could be saturated with 
available markers, and so that regions, referred to as 'coverage gaps/ could be 
pinpointed where markers were insufficiently informative in our sample to detect 
evidence of linkage. The Stage II screen followed up on regions flanking each 
marker that yielded peak lod scores approximately equal to or greater than the 
thresholds used for the coverage calculations, which were deemed regions of 
interest, and filled in coverage gaps. The results of the complete genome screen 
(Stages I and II) using 473 markers is described below. 

In addition, linkage disequilibrium analysis of an independently collected 
sample of 48 unrelated BP-I patients was conducted. These patients were from the. 
same ancestral population as the patients in the CR001 and CR004 pedigrees. LD 
analysis of further BP-I patients from the CRCV with markers in this 18q23 region 
is conducted to confirm and fine map a BP-I gene in this region. This approach, 
using additional BP-I patients from this CRCV population and additional markers 
identifies the region of maximum LD and can precisely localize a BP-I 
susceptibility gene. 

A conservative approach to linkage analysis was used in that almost all of 
the information for linkage is derived from individuals with a severe, narrowly 
defined phenotype. While this approach made it very unlikely that lod scores 
greater than conventional thresholds of statistical significance (e.g. _> 3) would 
be obtained, it provided confidence in the robustness of the most suggestive 
findings. 
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Three lines of evidence support the localization of a BP-I susceptibility 
locus to I8q22-q23: association analyses, linkage analyses, and direct observation 
of a conserved marker haplotype. 

Screening Lod Score Analyses 

Linkage analysis was performed using a narrow definition of the affected 
phenotype and a conservative model of genetic transmission. In particular, a high 
rate of phenocopies (non-genetic cases) was estimated. In an initial genome 
screen, two-point lod score analysis was conducted for all markers (Lathrop et al., 
1984; Ott, 1991; Terwilliger and Ott, 1994) and it was observed that all of the 
markers tested in 18q22-q23 displayed positive maximum iod scores 
(D18S64: 1.89 at recombination fraction (9)=0.18, D18S55: 1.45 at 9=0.18, 
D18S61: 1.75 at 9=0.16, and D18S70: 0.76 at 9 = 0.20). This suggestive 
evidence for linkage was obtained over a greater chromosomal length (at least 40 
centiMorgans (cM)) than in any other genome region. Eighteen additional 
microsatellite loci that cover the 18q22-q23 region at intervals of 0-5 cM were 
typed and linkage analyses were performed using the model employed for the 
initial genome screening studies. Fourteen of the markers, spaced across the 
region, displayed positive maximum lod scores (with six markers having lod 
scores > 1.0), with no peak localization. These results supported a localization of 
a BP-I gene in this region of chromosome 18. 

Association Analyses 

To further evaluate the evidence for this localization, an independent 
method was utilized which does not rely on specifying the mode of genetic 
transmission for BP-I. In isolated populations, it has previously been shown that 
the identification of genome regions where marker allele frequencies differ 
between patients and the background population can be used to map disease loci 
(Friedman et al., 1995). When such deviation is based on a substantially 
increased frequency of one or a few alleles in the patients, the region is almost 
certainly inherited identical by descent ("IBD"), with the disease gene, from 
common ancestors. The genome screening experiments followed the procedures 
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described by Boehnke, 1991 for directly estimating allele frequencies using 
genotypes from the family members in the study. For many of the markers tested 
in 18q23, the alleles most commonly observed in the BP-I patients are rare in the 
reference families of the Centre d'Etude du Polymorphisme Humain (CEPH), 
5 which are drawn from a variety of Caucasian populations (Dausset et al., 1990). 

For example, an allele at D18S70 of 124 basepairs in length was observed in 16 
out of 24 BP-I individuals tested but in only 3% of the chromosomes from the 
CEPH population. 

To evaluate whether these differences in allele frequency could simply be 
10 explained by genetic drift, a set of individuals sampled from the general Costa 

Rican population was genotyped using the 18q23 markers as well as a series of 
microsatellites (from other genome regions) that have been employed in 
comparative studies of several populations (Di Rienzo et al., 1994; Garza et al., 
1995). The Costa Rican reference sample showed no significant deviations from 
15 Hardy-Weinberg equilibrium for these markers, nor were there significant 

differences in allele frequency with other human populations (E. Rojas et al., 
unpublished observations). The marker allele frequencies for the patient sample 
were compared with the Costa Rican reference sample, without making any 
assumptions concerning linkages. For several of the 18q23 markers (D18S469, 
20 D18S554, D18S461, and D18S70), allele frequencies in the BP-I pedigree sample 

were highly significantly different from those in the Costa Rican population, but 
no such differences were observed for the markers from other genome regions. 

Revised Lod Score Calculations 

25 As allelic associations may profoundly affect the results of lod score 

calculations new linkage analyses using only affected individuals were carried out 
for the 18q23 markers, correcting for the observed significant associai is between 
BP-I and several of the marker alleles. The lod scores obtained for several of 
these markers (D18S380, D18S554, D18S461, and D18S70) provide independent 

30 support for the localization of a BP-I gene to this region, as indicated, for 

example, by marker D18S380 which yielded evidence of linkage but did not 
display significant allelic association with BP-I. 
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Conserved Marker Haplotypes 

The suggestion that a gene for BP-I is located in the 18q22-q23 region is 
further strengthened by the observation that BP-I segregates with particular marker 
haplotypes in both families (FIG. 1 and 5A and 5B). 18q22-q23 marker 
5 haplotypes for all genotyped individuals (and for deceased individuals whose 

genotypes could be reconstructed), BP-I segregates with particular marker 
haplotypes in both families (FIG. 1, 5A and 5B). In CR004, 16 of 17 affected 
individuals share portions of a marker haplotype from D18S64 to D18S70, a 
distance of about 40 cM (FIG. 5 A and 5B), with the majority of these individuals 

10 sharing at least 30 cM in common. Most affected individuals in CR004 share the 

distal segment of this haplotype, consisting largely of marker alleles that are rare 
in the general population of Costa Rica, indicating that the segment is inherited 
IBD from a common ancestor. A distinct, similarly sized haplotype is shared by 
seven of the nine BP-I individuals in CR001 (FIG. 5A). An identical haplotype, 

15 however, is observed in both families, in the region between D18S469 and 

D18S554 (a distance of about 6-7 cM). Of the BP-I individuals whose genotypes 
could be fully reconstructed, 22 out of 26 apparently share portions of this high- 
risk haplotype. Of the four BP-I individuals who do not display this haplotype 
sharing (FIG. 1, individuals IV-3, V-l, V-2. and V-12), three do not descend from 

20 a founding couple who lived in the 18th Century and were ancestral to the main 

branches of CR001 and CR004. As the interval between D185469 and D185554 
demonstrates sharing of a marker haplotype among most of the distantly as well as 
closely related BP-I individuals, the locus for BP-I susceptibility is probably within 
this six to seven cM segment. Probable recombinations in this interval in two 

25 individuals in CR004 (V-l 6 and V-l 7) were also identified. 

It is evident that the high risk haplotypes were introduced into the BP-I 
families through several ancestors, who were themselves distantly related to each 
other (FIG. 1). Although the BP-I patients studied are members of families that 
share common ancestry, the two pedigrees have been separated for at least seven 

30 generations. It is to be expected that such distantly related patients would share 

IBD a region of several cM surrounding a BP-I susceptibility gene but extremely 
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unlikely that they would share a segment of this size in any other area of the 
genome. 

Direct visual inspection is the best means of evaluating the evidence from 
multiple transmissions of a disease allele. Identifying such marker haplotypes 
facilitates mapping of disease loci when locus heterogeneity exists within extended 
pedigrees. The probability of locus (and ancestral) homogeneity increases when 
one samples families from relatively isolated populations. It is likely, however, 
that over several generations, many distantly related individuals will marry into 
such pedigrees, and thus the same disease allele (and a conserved marker 
haplotype surrounding the disease locus) may be introduced several times, and 
some individuals will be homozygous for rare alleles. 

Examination of the pedigrees studied clearly indicates multiple entries of a 
high risk haplotype in 18q22-q23. Portions of this haplotype occur in almost all 
patients in both pedigrees. Linkage analysis does not recognize the information 
conveyed by haplotypes shared via distant familial connections unless all relations 
are exactly specified and all markers are included in a single analysis; this would 
lead to unacceptable calculation times. In addition, such analyses generally do not 
sufficiently localize disease genes to permit positional cloning efforts. In contrast, 
identifying the region of maximal sharing of a marker haplotype in an isolated 
population can pinpoint the location of a gene for a complex disorder, even when 
the haplotype is observed in only a subset of patients (Nystrom-Lahti, et al., 
1994). Extended haplotypes, such as those observed, thus yield largely 
independent evidence for disease gene localization, corroborating the results of lod 
score and association analyses, and indicate regions to be highlighted in further 
mapping studies. 

Because the Costa Rican population is descended from only a small number 
of ancestors (M. Escamilla et al.. unpublished data) the localization of the BP-I 
susceptibility gene will be narrowed with additional BP-I samples from this 
population using linkage disequilibrium (LD) methods. Such an approach was 
recently used to follow up initial localization of a gene for insulin dependent 
diabetes mellitus (Davies et al., 1994; Copeman et al., 1995). For this reason 
examination of a sample of "unrelated" BP-I patients from the Costa Rican 
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population is expected to more precisely localize the putative BP-I gene by 
identifying the region of maximum LD (Lander and Shork, 1994; Puffenberger 
et al. 1994b; Collins, 1995). 

As predicted by prior genetic epidemiological studies, and by the 
segregation patterns of BP-I within extended pedigrees in which many apparent 
carriers are themselves unaffected, the high risk haplotype observed in CR001 and 
CR004 demonstrates incomplete penetrance (for example individual IV-26 is not 
affected with BP-I but apparently transmits illness to two sons). Although final 
diagnoses have been assigned to all identified members of CR001 and CR004 with 
BP-I and SAD-M, final diagnoses have not been to all members of these families. 
However, of those individuals who have been assigned final diagnoses other than 
BP-I and SAD-M, the shared marker haplotype is observed in all six of the 
individuals with BP-II (bipolar disorder with hypomania and major depression), 
MDD (major depressive disorder), or organic mood disorder. Six of eleven 
individuals with no apparent psychiatric diagnosis share this haplotype, however 
four of them transmit the haplotype (and presumably the risk of BP-I) to affected 
offspring. Precise delineation of penetrance in these families may not be possible 
until causative mutations are identified, as was recently demonstrated in pedigrees 
loaded for Hirschsprung's disease (Puffenberger et al. 1994b). 

The invention now being generally described, the following examples are 
provided for purposes of illustration only and are not to be considered to limit the 
invention. 

EXAMPLES 

Pedigrees 

Two independently ascertained Costa Rican pedigrees (CR001 and CR004) 
were chosen because they contained a high density of individuals with BP-I and 
because their ancestry could be traced to the founding population of the Central 
Valley of Costa Rica. The current population of the Central Valley (consisting of 
about two million people) is predominantly descended from a small number of 
Spanish and Amerindian founders in the 16th and 17th centuries Escamilla, M.A., 
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et al., (1996) Neuropsychiat. Genet. 67, 244-253. Studies of several inherited 
diseases have confirmed the genetic isolation of this population (Leon, P., et al. 
(1992) Proc. Natl, Acad. Sci. USA. 89,5181-5184, Uhrhammer, N., et al. (1992) 
Am. J. Hum. Genet. 57, 103-111). An extensive description of pedigrees CR001 
and CR004 is provided in a separate paper (Freimer, N.B., et al. (1996) 
Neuropsychiat. Genet. 67, 254-263). In the course of the study, two links 
between these pedigrees were discovered separately, however, because these links 
were discovered after the simulation analyses were completed and after the 
genome screening study had been initiated. 

All available adult members of these families were interviewed in Spanish 
using the Schedule for Affective Disorders and Schizophrenia Lifetime version 
(SADS-L) (Endicott, J. et al, (1978) Arch. Gen. Psych. 35, 837-844). Individuals 
who received a psychiatric diagnosis were interviewed again in Spanish by a 
research psychiatrist using the Diagnostic Interview for Genetic Studies (DIGS) 
(Nurnberger, J.L. et al. (1994) Arch. Gen. Psychiat. 51, 849-859). This recently 
developed instrument is similar to, but more detailed than SADS-L. The 
interviews and medical records were then reviewed by two blinded best estimators 
who reached a consensus diagnosis. The diagnostic procedures are described in 
detail in Freimer, N.B., et al. (1996) Neuropsychiat. Genet. 67, 254-263 
(incorporated by reference herein). 

Unrelated CRCV BP-I Patient Study 

BP localizations obtained through the CRCV pedigree studies were 
confirmed by genotyping an independently collected sample of 48 unrelated BP-I 
patients from the CRCV. In this fine mapping LD analysis, 48 unrelated BP-I 
patients from the CRCV were identified and genotyped using microsatellite 
markers spaced at narrow intervals across chromosome 18. As these patients are 
descended from the same ancestral population as the patients in the pedigrees 
previously studied (CR001 and CR004), many of them should share disease 
susceptibility alleles inherited identical by descent (IBD) from one or a few 
common ancestors, and linkage disequilibrium (LD) should be present at marker 
loci surrounding the disease genes. 
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The sample of 48 BP-I patients included 25 women and 23 men who were 
recruited from psychiatric hospitals and clinics in the CVCR. These patients were 
ascertained only on the basis of diagnosis and CV ancestry, and were not selected 
on the basis of history of BP illness in family members. A structured interview of 
each patient was conducted by a psychiatrist, and medical and hospital records 
were collected. Ascertainment and diagnostic procedures were as described above. 
However, in order to further lessen the probability of phenocopies among this 
unrelated sample, for which we lacked pedigree information, the affected 
phenotype was defined even more narrowly than in the pedigree study. Individuals 
considered affected in this study had to have suffered at least two disabling 
episodes of mania (requiring hospitalization) and a first onset of the illness before 
age 45. 

Genealogical research on each of the 48 BP-I patients confirmed that on 
average. 70% of their great-grandparents were born in the CRCV. Individuals 
whose great-grandparents were bom in the CRCV were considered likely to be 
descended from the original Spanish and Amerindian founders of the CRCV. 
Genealogical research showed that 2 patients are first cousins and the remaining 46 
have no relationship within the past 4 generations. 

Genotyping Pedigree Studies 

Linkage simulations were used to select the most informative individuals 
from pedigrees CR001 and CR004 for genotyping studies (Freimer. N.B., et al. 
(1996) Neuropsychiat. Genet. 67, 254-263). Under a 90% dominant model, 
simulation analyses with these individuals suggested that evidence of linkage would 
likely be detected (e.g. a probability of 92% of obtaining lod > 1.0 in the 
combined data set) using markers with an average heterozygosity of 0.75 spaced at 
10 cM intervals (as discussed in 24)). For the Stage I screen, the most 
polymorphic markers (307 in total) were chosen, placed at approximately 10 cM 
intervals on the 1992 Genethon map (Houwen, R., et al. (1992) Nature 359, 794- 
801). These markers were then supplemented by a small number of markers from 
the Cooperative Human Linkage Center (CHLC) public database. For the Stage II 
screen. 166 markers were added from newer Genethon and CHLC maps as they 
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became available (Murray. J.C. et al. (1994) Science 265, 2049-2054, Gyapay, 
G., et al. (1994) Nature Genet. 7,246-339) and from the public database of the 
Utah Center for Genome Research, all of which are publicly available. DNA 
samples (from individuals in the CEPH families) that were used for size standards 
for Genethon and CHLC markers were included in the experiments to permit 
comparison of allele sizes between members of the CRCV population and 
individuals in the CEPH database. Marker D18S1009 is a previously unpublished 
Genethon marker, however information about this marker is now available from 
the Genome Database. 

Genotyping procedures were as described previously DiRienzo, A. et al. 
(1994) Proc. Natl. Acad. Sci. USA 91, 3166-3170 (incorporated by reference 
herein). Briefly, one of the two PCR primers was labeled radioactively using a 
polynucleotide kinase and PCR products were run on polyacrylamide gels. 
Autoradiographs were scored independently by two raters. Data for each marker 
were entered into the computer database twice and the resultant files were 
compared for discrepancies. 

Association Analysis 

Associations between markers on 18q23 and bipolar disease were assessed 
by comparing the frequency of marker alleles as estimated from affected 
individuals in the families with the frequency of marker alleles in the Costa Rican 
reference population. Marker allele frequencies were first estimated on a 
combined data set including the bipolar patients in the pedigrees (with correction 
for dependency due to family relationships as per Boehnke 1991) and the 
population sample. Alleles that were rare in the combined data set were collapsed 
into a single allele for each marker, until no alleles with a frequency of less than 
4% remained, to prevent chance occurrence of rare alleles from having 
disproportionate effects in subsequent comparisons. Following this recoding, 
allele frequencies were again estimated on the combined data. Marker allele 
frequencies were then estimated separately in the reference population and in 
affected family members. A similar procedure was used by Schellenberg et al. 
(1987). The likelihood obtained in analyzing the combined data was compared 
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with the product of the likelihoods obtained in the separate analyses, to test the 
null hypothesis of no difference between allele frequencies in the Costa Rican 
reference sample and affected individuals from the BP-I families. Under the null 
hypothesis, this likelihood ratio statistic is distributed as a Chi-square random 
variable with n-l degrees of freedom, where n indicates the number of alleles for 
the marker. 

Haplotype Construction 

Haplotypes for the pedigrees in 18q were constructed by hand, without 
knowledge about the diagnostic status of family members. A minimum 
recombination strategy was the parsimonious method chosen to guide haplotype 
construction. Three individuals who carry diagnoses of BP-I or SAD-M are 
deceased but their haplotypes could be partially (IV-12) or nearly fully 
reconstructed (IV-9, IV- 18). The haplotype information is presented schematically 
in FIG. 5A and 5B. Only BP-I and SAD-M individuals and their direct ancestors 
are depicted. As indicated in the discussion of the conserved haplotypes, there is 
no adequate means to statistically represent the evidence deriving from such 
transmission. 

Definition Of The Boundaries Of The Candidate Region By Additional 
Marker Typing Studies 

Statistical analysis suggested a BP-I localization in 18q23-qter, but the 
strongest evidence derives from an extended marker haplotype that is shared by 
most of the BP-I patients in the pedigree study, particularly in Family CR004. 
This haplotype also helps define a candidate region for the additional mapping. 

Genetic approaches, as described below, are used to aim for a BP-I 
candidate region of less than one Mb. As new markers are tested, and increased 
haplotype sharing is identified, one moves closer to the causative gene 
(Puffenburger, 1994a). 

The most likely region for a BP-I locus is between markers D18S469 and 
D18S554, a distance of 7 cM. This location is suggested by recombination events 
that appear to interrupt the conserved haplotype in Family CR004. Most patients 
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share a haploiype from D18S554 towards the telomere (marker D18S70), and the 
majority share a hapiotype from above D18S469 to D18S70. However, two 
patients share only the centromeric portion of the hapiotype (a total of 12 markers 
covering more than 15 cM); that is, there is a "break" below D18S469. 
5 Additional possible evidence for this location is provided by two of the three 

affected individuals in CR0O4 who do not share the extended hapiotype observed 
in the other patients. These individuals may share this hapiotype at markers 
D18S469 and D18S554. 

An additional six markers from Genethon have been mapped within 
10 D18S469 and D18S554 dividing the interval into one to two cM segments. A 

further set of new markers has also been mapped to this region, by the CHLC 
(1994). All of the available Genethon and CHLC markers are used to genotype 
DNA from the previously investigated individuals from CR001 and CR004. From 
these experiments a refined region of maximal hapiotype sharing is delineated. In 
15 addition, IBD and IBS sharing is distinguished for the patients who currently 

possess "BP-I" alleles only for markers D18S469 and D18S554. 

It is expected that the above described experiments will delineate a likely 
candidate region of as small as one cM, permitting initiation of physical mapping 
and cloning studies. However, it is necessary to enlarge the study sample and to 
20 identify additional markers in order to proceed with the detailed mapping of this 

region. 

For all of the genotyping studies used to characterize a putative BP-I locus 
in the 18q23-qter chromosomal region, standard procedures (Di Rienzo et al., 
1994) are used. One of the PCR primers is radioactively labeled with P32 using 
25 T4 kinase, and the PCR products are size separated using denaturing sequencing 

gels and detected by autoradiography. The sizes of alleles are determined to the 
level of a single base pair by comparison with known standards (Di Rienzo et al., 
1994). 

30 
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Genotyping of Unrelated BP-I CRCV Patients 
Twenty-seven markers were used to genotype all 48 individuals (as well as 
53 individuals used to establish genetic phase) at approximately 5 cM intervals 
along the entire chromosome 18. It was hypothesized that such a screen would 
5 permit the evaluation of evidence in the 18q22-q23 region and also to investigate 

other regions on chromosome 18 in which linkage to BP has been suggested by 
other groups in other populations. For each individual, two-marker haplotypes in 
each of the 26 inter- marker intervals were investigated. For 38 of the 48 BP-I 
patients, genotypes of parents or children were available to assist in phase 

10 determination. Because of phase ambiguities in the remaining 10 individuals, 

minimal and maximal two-marker haplotype sharing was evaluated as follows: (1) 
Minimal: the number of individuals (and chromosomes) who definitely shared a 
chromosomal segment defined by a particular pair of alleles (phase known 
chromosomes) and (2) Maximal: the number of individuals (and chromosomes) 

15 who could possibly share a chromosomal segment defined by a particular pair of 

alleles (includes phase unknown chromosomes). The threshold used to identify 
areas of high IBD sharing of chromosomes in this initial screen was designated as 
maximal sharing of a two-marker Kaplotype by 50% or more of the 48 individuals 
(or 25% or more of the 96 chromosomes). 

20 Arbitrary thresholds were designated to identify possible areas of high IBD 

sharing among the 48 patients. Eight of the 26 regions passed this screen. Within 
each of these 3 regions, one to three additional markers were typed to permit 
detection of LD, if present, over regions of one to two cM. 

A total of 42 chromosome 18 markers were used to genotype the study 

25 sample: 

D18S59, D18S1140, D18S476, D18S481. D18S391, D18S452, D18S843, 
D18S464, D18S1153, D18S378, D18S53. D18S453, D18S40, D18S66, D18S56, 
D18S57, D18S467, D18S460, D18S450, D18S474, D18S69, D18S64, D18S1134, 
D18S1147, D18S60, D18S68, D18S55, D18S477, D18S61, D18S488, D18S485, 

30 D18S541, D18S870, D18S469, D18S874, D18S1121, D18S380. D18S1009, 

D18S844, D18S554, D18S46L D18S70 (from pter to qter). Of these 42 markers, 
seven are located within the 18q22-23 region extending from the marker D18S469 
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to marker D18S554 (inclusive). This region is referred to as the 18q22-23 region. 

For each marker the likelihood that a particular allele (or alleles) is over- 
represented on disease chromosomes, as compared to non-disease chromosomes 
was evaluated. The results of this likelihood test provide a conservative but 
powerful measure of LD between two loci. 



Pedigree Statistical Analyses 

Two-point linkage analyses were performed for all markers. Marker allele 
frequencies were estimated from the combined data set with correction for 
dependency due to family relationships (Boehnke, M. (1991) Am. J. Hum. Genet. 
48, 22-25). The linkage analyses for Stages I and II included the 65 individuals 
who were genotyped as well as an additional 65 individuals who had been 
diagnostically evaluated but not genotyped. Only individuals with BP-I were 
considered affected with the exception of two persons, one in each family, who 
carry diagnoses of schizoaffective disorder manic type (SAD-M), The SAD-M 
individuals were included as affected because BP-I and SAD-M are often difficult 
to distinguish from each other based on their clinical presentation and course of 
illness (Goodwin, F.K. et al. (1990) in Manic Depressive Illness. (Oxford 
University Press, New York), pp. 373-401; Freimer, N.B et al. (1993) in The 
Molecular and Genetic Basis of Neurological Disease, pp. 951-965; Freimer, N.B. 
et al. (1996) Neuropsychiat. Genet. 67, 254-263; and Freimer, N.B. et al (1996) 
Nature Genetics 12:436-441, all incorporated by reference herein). In all, 
20 individuals were designated as affected within CR004 (Copeman, J.B., et al. 
(1995) Nature Genet. 9, 80-85 available for genotyping) and 10 individuals from 
CR001 (Kelsoe, J.R. et al. (1989) Nature 342, 238-243 available for genotyping). 
The phenotype for all other individuals was designated as unknown except for 17 
individuals who were designated as unaffected because they had been thoroughly 
clinically evaluated, showed no evidence of any psychiatric disorder, and were 
well beyond the age of risk (50) for BP-I (linkage simulation studies indicated that 
these unaffected individuals contributed little information to the linkage analysis). 
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Linkage analyses were performed using a nearly dominant model (assuming 
penetrance of 0.81 for heterozygous individuals of 0.9 for homozygotes with the 
disease mutation). This model was chosen from five different single-locus models 
(ranging from recessive to nearly dominant) due to its consistency with the 
segregation patterns of BP in the two pedigrees and because it had demonstrated 
the greatest power to detect linkage in simulation studies Freimer, N.B., et al. 
(1996) Neuropsychiat. Genet. 67, 254-263. Based on Costa Rican epidemiological 
surveys Escamilla, M.A., et al., (1996) Neuropsychiat. Genet. 67, 244-253, the 
population prevalence of BP-I was assumed to be 0.015 (and thus the frequency of 
the disease allele was assumed to be 0.003)(based on epidemiological surveys in 
Costa Rica, Adis, G. (1992) "Disordenes mentales en Costa Rica: Observaciones 
Epidemiologicas," (San Jose, Costa Rica: Editorial Nacional de Salud y Seguridad 
Social)). The frequency of BP-I in individuals without the disease allele was 
conservatively set at 0.01 which effectively specified a population phenocopy rate 
15 of 0.67 (i.e. an affected individual in the general population has a 2/3 probability 

of being a phenocopy). For multiply affected families, the probability that a gene 
is segregating is highly increased, which implies that affected individuals in our 
study pedigree have a lower probability to be phenocopies than affected individuals 
in the general population, particularly those with several affected close relatives 
20 (the exact probabilities are dependent on the degree of relationship between patents 

and the number of intervening unaffected individuals). These parameters were 
chosen to ensure that most of the linkage information derives from affected 
individuals. The rationale for selecting these parameters and results of analyses 
that demonstrate the conservatism of this model are described by Freimer, N.B., 
et al. (1996) Neuropsychiat. Genet. 67, 254-263. The LINKAGE package 
(Lathrop et al., (1984) Proc. Natl. Acad. Sci. USA 81. 3443-3446) was used for 
lod score analysis and to obtain maximum likelihood estimates of the marker allele 
frequencies, taking into account the existing family relationships (see Boehnke, 
Am. J. Hum. Gent. 48, 22-25 (1991)). 



25 
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A likelihood test of disequilibrium (J. Terwilliger, Am. J. Hum. Genet. 56, 
777 (1995)) was used to estimate a single parameter, lambda, that quantifies the 
over-representation of marker alleles on disease chromosomes as compared to non- 
disease chromosomes. We chose this method of analysis over another commonly 
used disequilibrium analysis method, the transmission disequilibrium test (TDT, R. 
Spielman et al.. Am. J. Hum. Genet. 52, 506 (1993)) because data from all 48 
BP-I patients could be used in the likelihood approach. Effective use of the TDT 
requires phase-known, heterozygous parental chromosomes. We do not have 
parental genotypes for 20 of the 48 BP-I patients. Simulations indicated that with 
our data, the likelihood test of disequilibrium would be more powerful than the 
TDT. Lambda has been shown to be a superior measure for LD fine mapping, 
compared to other frequently used measures, because it is directly related to the 
recombination fraction between the disease and the marker loci. Non-disease 
chromosomes were chosen from the phase-known chromosomes of parents, 
spouses and children of affected individuals, if available. Designation of 
chromosomes of family members as non-disease in a disorder such as BP-I, which 
is not fully penetrant, necessitates specifying a model of disease transmission. The 
same model of transmission was employed in this LD likelihood test as was used 
in the initial genome screen of the pedigrees CR001 and CR002 described herein. 
One parameter was specified differently from the genome screen: the phenocopy 
rate was set to zero in the LD likelihood analysis. A phenocopy rate was not 
specified in the transmission model because the effect of phenocopies will be 
"absorbed' by the lambda parameter, in that presence of phenocopies in our 
sample will serve to erode the association between marker alleles and disease, and 
hence reduce the estimate of lambda. 

Coverage 

To access coverage for a marker, the number of informative meioses at the 
estimated recombination fraction was calculated using the estimate of the variance 
(the inverse of the information matrix) (Petrukhin, K.E. et al. (1993) Genomics 
15, 76-85). Alternatively, when the estimated frequency of recombination was 
close to 0 or 1, Edwards' equation was applied to calculate the equivalent number 
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of observations (Edwards. J.H. (1971) Ann. Hum. Genet. 34, 229-250). These 
meioses represent the amount of linkage information provided by the marker, 
given the pedigree structure and the genetic model applied. Linkage to the marker 
in question was then assumed and the load score that would be observed as a 
disease gene is hypothetically moved in increments away from that marker was 
calculated. All regions abound a marker that would have generated a lod score 
that exceeded our thresholds for possible linkage (0.8 in CR001, 1.2 in CR004, 
and 1.6 in the combined data) were considered covered. These lod score 
thresholds were derived from simulation analyses showing the expected 
distribution of lod scores under linkage and non-linkage Freimer, N.B., et al. 
(1996) Neuropsychiat. Genet. 67, 254-263, and approximately represent a result 
that is 250 times more likely to occur in linked simulations than in unlinked 
simulations. Coverage maps were constructed (FIG. 5 A and 5B) by 
superimposing the regions covered by each market on the genetic map of each 
15 chromosome. At the end of the Stage II screen, a total of 473 microsatellite 

markers had been typed with genome coverage (in the combined data set) of over 
94%. Possible coverage gaps are indicated by unshaded areas and are mainly 
concentrated near telomeres. Because the coverage calculations make use of 
market informativeness within the pedigrees, the coverage approach thus permits 
20 detection of instances where markers with expected high heterozygosities are 

uninformative in our data set. 

Pedigree Linkage Analysis Results 

Of the 473 microsatellites analyzed with two-point linkage tests, 23 
25 markers exceeded the empirically determined thresholds designated for the 

coverage calculations (in either CR001, CR004, or in the combined data set). The 
location of these markers, the peak lod scores obtained in each family and in the 
combined data set, and the maximum likelihood estimate of the recombination 
fraction (0) at which these lod scores were observed are indicated in Table 1. The 
30 approximate chromosomal locations of these markers are also depicted in FIG. 5A 

and 6B. The distribution of lod scores (for the maximum likelihood estimate of 0 
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in the combined data set) across the genome is displayed by chromosome in FIG. 
2. 

Five markers in the 18q22-q23 region exceeded the lod score thresholds in 
pedigree CR004. See FIG. 2 and 3. 

5 

Unrelated BP-I CRCV Patient Study Results 
Out of the forty-two markers tested, eight displayed evidence of over- 
representation of a particular allele on disease chromosomes. Eight of the 42 
markers had -2*ln(likelihood ratio) statistics > 1,0. Three other markers had - 
10 2*ln(likelihood ratio) statistics >0 and <0.62. The results are shown in Table I: 

Table I 



Marker 


Allele Size 


Frequency on 
non-disease 
Chromosomes 


Frequency on 

Disease 
Chromosomes 


D18S59 


154 


0.121 


0.572 


D18S476 


271 


0.470 


0.771 


D18S467 


172 


0.384 


0.693 


D18S61 


177 


0.074 


• 0.326 


D18S485 


182 


0.237 


0.586 


D18S870 


179 


0.405 


0.657 


D18S469 


234 


0.128 


0.450 


D18S1121 


168 


0.171 


0.553 



25 Five of these eight marker loci were clustered within a small region of 18q22-q23 

that overlapped the region of maximal haplotype sharing among affected 
individuals in CR001 and CR004, The strongest LD in 18q22-q23 was observed 
at D18S1121 (-2*In(likelihood ration) of 5.03, p=0.01); it is noteworthy that at 
this locus a specific allele (of 168 bp) is over- represented on the disease 

30 chromosomes in the population BP-I sample and is also present on the predominant 

high-risk haplotype within the pedigrees. Fifteen of 17 affected individuals in 
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CR004, and 4 of 9 affected individuals in CR001 possessed the 168 b.p. allele at 
D18S1121. 

The LD and pedigree findings in the 18q23 region denote a clearly 
delineated region flanked by markers D18S469 and D18S554 that contains a BP-I 
susceptibility locus. Results of additional LD analyses using additional markers 
are suggestive of a potential gene involved in BP-I susceptibility in a region of 
18q22-q23 defined by markers D18S1121 and D18S380, inclusive. These regions 
are distinct from other regions on chromosome 18 that have been suggested as 
linked to mood disorder phenotypes (more broadly defined than BP-I). See FIG. 
5A, 5B. In contrast to previous reports by Berrettini et al. and Stine et al., 
suggesting possible linkage between mood disorder and markers in the 
pericentromeric region of chromosome 18, our results did not show any evidence 
for association of BP-I with any pericentromeric markers (D18S378, D18S53, 
D18S453 or D18S40). 



Additional Fine Mapping Studies 

In order to further delineate regions of BP-I susceptibility within the 18q22- 
q23 region, additional unrelated BP-I patients from the CRCV as well as other 
populations are diagnosed and genotyped both with the markers described herein 

20 as well as additional markers in the 18q23 region that are known as well those yet 

to be identified. Additional markers are available from the Cooperative Human 
Linkage Center (CHLC) public database, from newer Genethon and CHLC maps 
as they become available (Murray, J.C. et al. (1994) Science 265, 2049-2054, 
Gyapay, G., et al. (1994) Nature Genet. 7,246-339) and from the public database 

25 of the Utah Center for Genome Research (all of which are incorporated by 

reference herein). The web addresses for Genethon and CHLC are: Genethon 
(http://www. genethon. fr/genethon_en.html), CHLC 

(http://gopher.chlc.org/HomePage.html). These databases are ail linked, and one 
of ordinary skill in the art can readily access the information available from these 
30 databases. 

The markers shown in FIG. 6B. can be used to genotype the CRCV 
pedigrees and unrelated BP-I patients described herein as well as other BP-I 
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affected individuals and pedigrees. See FIG. 6B (portion of a chromosome 18 
map available from the Whitehead Institute, web address: 
http://133. 30.8. 1: 8080/ = ® = : www-genome. wi.mit.edu. (incorporated herein by 
reference)). The fine mapping techniques described herein in conjunction with the 
teachings regarding the 18q23 region can be used to further narrow the BP-I 
susceptibility region. 



Identification Of New Markers From the 18q23 Region 

Cloned human genomic DNA covering the target region is assembled. 
Microsatellite sequences from these clones are identified. A sufficient area around 
the repeat to enable development of a PCR assay for genomic DNA is sequenced, 
and it is confirmed that the microsatellite sequence is polymorphic, as several 
uninformative microsatellites are expected in any set. Several methods have been 
routinely used to identify microsatellites from cloned DNA, and at this time no 
single one is clearly preferable (Weber, 1990, Hudson et a!., 1992). Most of 
these require screening an excessive number of small insert clones or performing 
extensive subcloning using clones with larger inserts. 

New strategies have recently been developed which permit the use of the 
several different microsatellites to be found within a single large insert clone 
without requiring extensive subcloning. A method for direct identification of 
microsatellites from yeast artificial chromosomes (YACs) provides several new 
markers from the target region. This procedure is based on a subtractive 
hybridization step that permits separation of the target DNA from the vector 
background. This step is useful because the human DNA (the YAC) constitutes 
only a small proportion of the total yeast genomic DNA. 

YAC clones (with inserts averaging about 750 Kb of human genomic DNA) 
that span the 18q22-q23 region have already been identified by the 
CEPH/Genethon consortium (Cohen et al., 1993) and are publicly available. The 
markers from YACs that have been mapped to portions of the candidate region 
that are not well represented by currently available markers are first isolated. By 
.typing these markers in the families and the M LD" sample, as described above, it 
is possible to narrow the candidate region, perhaps to a size of less than one to 
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two cM, thus permitting limitation of the segment in which more extensive 
mapping efforts are applied. 

Briefly, the microsatellite identification procedure is performed as follows: 
A subtractive hybridization is performed using genomic DNA from a target YAC 
together with an equivalent amount of a control DNA. This procedure separates 
the YAC DNA from that of the yeast vector. Following the subtraction procedure 
the subtracted YAC DNA is purified, digested with restriction enzymes and cloned 
into a plasmid vector (Ostrander et al., 1992). The cloned products of each YAC 
are screened using a CA(15) oligonucleotide probe. Each positive clone (i.e. those 
that contain TG-repeats) is sequenced to identify primers for PCR to genotype the 
BP-I samples. 

An alternative approach, based on using a set of degenerate sequencing 
primers that anneal directly to the repeat sequence, permitting direct thermal cycle 
sequencing (Browne & Litt, 1992), can also be used. 

Once the candidate region is narrowed to a size of less than about 500 to 
1000 Kb. a contiguous array (contig) of clones with smaller inserts than YACs, 
mainly PI clones, is developed. PI clones are phage clones specially designed to 
accommodate inserts of up to 100 Kb (Shepherd et al., 1994). 

Development Of A Physical Map Of The 18q22-q23 Region 

In parallel with the genetic mapping, a physical map of the 18q22-q23 
region is developed. The backbone of this effort is the assembly of contigs of 
large insert clones. Low resolution contigs for most of the human genome are 
already available using the YACs developed by CEPH (Cohen et al., 1993). 
Although these have been individually verified and checked for overlap with other 
YACs, there is a high rate of chimerism in the YACs and insufficient evidence to 
definitively confirm the order of the YACs. In addition, because of their large 
size these YACs are particularly cumbersome to work with. Nevertheless, they 
provide a useful framework to start constructing high resolution contigs. 

Once a candidate region of less than about five cM is delineated, the 
studies to develop a physical map are commenced. Because of the disadvantages 
of relying solely on YACs, and because positional cloning is facilitated by the 
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availability of a higher resolution map, contigs are generated using PI clones once 
the candidate region is narrowed to less than one Mb, by LD mapping in the 
expanded population sample using the new markers identified from the YACs. 

Once a region of 500-1000 Kb or less is defined, physical mapping and 
cloning are computed using PI clones rather than YACs, and PI contigs over such 
a region are constructed. The Pis are used to identify additional markers for the 
further positional cloning steps as well as the screening for rearrangements. 

The starting point of contig construction is the microsatellite sequences and 
non-polymorphic STSs that derive from the few YACs that surround the 
genetically determined candidate region. These STSs are used to screen the PI 
library. The ends of the Pis are cloned using inverse PCR and used to order the 
Pis relative to each other. Amplification in a new PI will indicate that it overlaps 
with the previous one. Fluorescent in situ hybridization (FISH) permits ordering 
of the majority of the Pis (Pinkel, 1988; Lichter, 1991). The original set of Pis 
serves as building blocks of the complete contig; each end clone is used to re- 
screen the library and in this way Pis are added to the map. 

From each Pi additional microsatellites are identified as previously 
described. This allows further reduction of the candidate region. When the 
region is narrowed to less than one Mb in size, positional cloning efforts are 
initiated. 

Use Of PI Clones To Identify Candidate cDNAs For Screening 
For Mutations In The DNA Of BP-I Patients 

The Pi clones are used to identify candidate cDNAs. The candidate 
cDNAs are subsequently screened for mutations in DNA from BP-I patients. 
From the minimal candidate region defined by genetic mapping experiments a 
segment is left that is sufficiently large to contain multiple different genes. 

Identification Of Coding Sequences 

Coding sequences from the surrounding DNA are identified, and these 
sequences are screened until a probable candidate cDNA is found. Much of the 
human genome will be sequenced over the next few years, in which case it may 
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become feasible to identify coding sequences through database screening. 
Candidates may also be identified by scanning databases consisting of partially 
sequenced cDNAs (Adams et al., 1991), known as expressed sequence tags, or 
ESTs. These resources are already largely developed, and include upwards of 
5 100,000 cDNAs, the majority expressed primarily in the brain. It is not yet clear, 

however, that the complete set of cDNAs will be mapped to specific chromosomal 
locations in the near future, and that their data will soon be made publicly 
available. The database can be used to identify all cDNAs that map to the 
minimal candidate region for BP-I. These cDNAs are then used as probes to 

10 hybridize to the PI contig, and new microsatellites are isolated, which are used to 

genotype the "LD" sample. Maximal linkage disequilibrium in the vicinity of one 
or two cDNAs is identified. These cDNAs are the first ones used to screen 
patient DNA for mutations. Database screening has already been used to identify 
a gene responsible for familial colon cancer (Papadopolous et al., 1993). 

15 Coding sequences are also identified by exon amplification (Duyk et al., 

1990; Buckler et al., 1991). Exon amplification targets exons in genomic DNA by 
identifying the consensus splice sequences that flank exon-intron boundaries. 
Briefly, exons are trapped in the process of cloning genomic DNA (e.g. from Pis) 
into an expression vector (Zhang et al., 1994). These clones are transfected into 

20 COS cells, RT-PCR is performed on total or cytoplasmic RNA isolated from the 

COS cells using primers that are complementary to the splicing vector. Exon 
amplification is tedious but routine; for example, the system developed by Buckler 
et al. (1991). This method is probably preferable to another widely used 
approach, direct selection, which involves screening cDNAs using large insert 

25 clone contigs, with several steps to maximize the efficiency of hybridization and 

recovery of the appropriate hybrid (Lovett et al., 1991). Although direct selection 
is more efficient than exon amplification (Del Mastro et al., 1994), it may not be 
practical as it depends on the candidate cDNA being expressed in the tissue from 
which the cDNA library was made; there is no prior information to indicate the 

30 tissue or developmental stage in which BP-I genes would be expressed. 

Once cDNAs are identified the most plausible candidates are screened by 
direct sequencing. SSCP or using chemical cleavage assays (Cotton et al. 1988). 
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The data are also evaluated for clues to the possible identity or mode of 
action of BP-I mutations. For example, it is known that trinucleotide repeat 
expansion is associated with the phenomenon of anticipation, or the tendency for a 
phenotype to become more severe and display an earlier age of onset in the lower 
generations of a pedigree (Ballabio, 1993). Several investigators have suggested 
that segregation patterns of BP-I are consistent with anticipation (Mclnnis et al., 
1993; Ny lander et al., 1994). The apparent transmission of BP-I, in association 
with the conserved 18q23 haplotype is constant with anticipation. Therefore, once 
the candidate region is narrowed to its minimal extent, the PI clones are screened 
using trinucleotide repeat oligonucleotides (Hummerich et al., 1994). A PCR 
assay is developed and patient DNAs are screened for expanded alleles. 

Genetic and physical data help to map the bipolar mood disorder gene to 
the 18q22-q23 region of chromosome 18. New markers from this region are 
tested in order to locate the bipolar mood disorder gene in a region small enough 
to provide higher quality genetic tests for bipolar mood disorder, and to 
specifically find the mutated gene. Narrowing down the region in which the gene 
is located will lead to sequencing of the bipolar mood disorder gene as well as 
cloning thereof. Further genetic analysis employing, for example, new 
polymorphisms flanking D18S59 and D18S476 as well as the use of cosmids, yeast 
artificial chromosome (YAC) clones, or mixtures thereof, are employed in the 
narrowing down process. The next step in narrowing down the candidate region 
includes cloning of the chromosomal region of 18q22-q23 including proximal and 
distal markers in a contig formed by overlapping cosmids and YACS. Subsequent 
subcloning in cosmids, plasmids or phages will generate additional probes for 
more detailed mapping. 

The next step of cloning the gene involves exon trapping, screening of 
cDNA libraries, Northern blots or rt PCR (reverse transcriptase PCR) of samples 
from affected and unaffected individuals, direct sequencing of exons or testing 
exons by SSCP (single strand conformation polymorphism), RNase protection or 
chemical cleavage. 

Flanking markers on both sides of the bipolar mood disorder gene 
combined with D18S59 and D18S476 or a number of well-positioned markers that 
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cover the chromosomal region (18q22-q23) carrying the disease gene, can give a 
high probability of affected or non-affected chromosomes in the range of 80-90% 
accuracy, depending on the informativeness of the markers used and their distance 
from the disease gene. Using current markers linked to bipolar mood disorder, 
5 and assuming closer flanking markers will be identified, a genetic test for families 

with bipolar mood disorder will be for diagnosis in conjunction with clinical 
evaluation, screening of risk and carrier testing in healthy siblings. In the future, 
subsequent delineation of closely linked markers which may show strong 
disequilibrium with the disorder, or identification of the defective gene, could 
10 allow screening of the entire at-risk population to identify carriers, and provide 

improved treatments. 

Treatment of BP-I Patients Using Genotype Data 

Using the fine mapping techniques described herein, BP-I susceptibility loci 

15 or genes in the 18q22-q23 region are identified and used to genotype patients 

diagnosed phenotypically with BP-I. Genotyping with the markers described 
herein as well as additional markers permits confirmation of phenotypic BP-I 
diagnoses or assist with ambiguous clinical phenotypes which make it difficult to 
distinguish between BP-I and other possible psychiatric illnesses. A patient's 

20 genotype in the 18q22-q23 region is determined and compared with previously 

determined genotypes of other individuals previously diagnosed with BP-I. Once 
an individual is genotyped as having a BP-I susceptibility locus in the 18q22-q23 
region, the individual is treated with any of the known methods effective in 
treating at least certain individuals affected with BP-I. These known methods 

25 include the administration of drugs including antidepressants such as lithium salts, 

carbamazepine. valproic acid, lysergic acid diethylamide (LSD), p- 
chlorophenylalanine, p-propyldopacetamide dithiocarbamate derivatives such as 
FLA 63; anti-anxiety drugs, including diazepam, monoamine oxidase (MAO) 
inhibitors including iproniazid, cloryline, phenelzine and isocarboxazid; biogenic 

30 amine uptake blockers, e.g., tricyclic antidepressants such as desipramine, 

imipramine and amitriptyline; serotonin reuptake inhibitors, e.g., fluoxetine; 
antipsychotic drugs such as phenothiazaine derivatives (e.g., chiorpromazine 
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(thorazine), and trifluopromazine, butyrophenones (e.g., haloperidol (Haldol), 
thioxanthine derivatives (e.g., chlorprothixene); and dibenzodiazpines (e.g., 
chlozapine); benzodiazpines; dopaminergic agonist and antagonists, e.g., L-Dopa, 
cocaine, amphetamine, a-methyl-tyrosine, reserpine, tetrabenazine, benzotropine, 
pargyline; noradrenergic agonists and antagonists, e.g., clonidine, 
phenoxybenzamine, phentolamine, and tropolone. Many of these drugs are used 
in combination. 

Studies are conducted correlating effective treatments with BP-I genotypes 
in the 18q22-q23 region to determine the most effective treatments for particular 
genotypes. BP-I patients can then be genotyped in the 18q22-q23 region and the 
statistically most effective treatment can be determined as a first course of therapy. 

All publications and patent applications mentioned in this specification are 
herein incorporated by reference to the same extent as if each individual 
publication or patent application was specifically and individually indicated to be 
15 incorporated by reference. 

The invention now being fully described, it will be apparent to one of 
ordinary skill in the art that many changes and modifications can be made thereto 
without departing from the spirit or scope of the appended claims. 
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1. A method of predicting a patient's likelihood of developing bipolar mood 
disorder comprising 

obtaining a sample of said patient's DNA, 

determining said patient's genotype in a region on the long arm of 
chromosome 18, comprising determining said patient's allele sizes at markers 
located between markers D18S469 and D18S554, inclusive, and 

comparing said patient's genotype to genotypes of affected individuals. 

2. The method of claim 1, wherein said patient's genotype is determined in a 
region on the long arm of chromosome 18 between markers D18S1121 and 
D18S380 t inclusive. 

3. The method of claim 1, wherein said markers consist essentially of markers 
located between D18S1121 and D18S380. 

4. The method of claim 1, wherein said genotypes of affected individuals 
comprise allele sizes depicted in Fig. 5A and 5B. 

5. The method of claim 1, wherein said patient has Spanish or Amerindian 
ancestry. 

6. A method of predicting a patient's responsivity to drug treatment for 
bipolar mood disorder comprising: 

determining said patient's genotype in a region on the long arm of 
chromosome 18, wherein said genotype comprises said patient's allele sizes at 
markers located between markers D18S469 and D18S554, inclusive, 

comparing said patient's genotype to genotypes of affected individuals, 
wherein said genotype comprises said affected individuals' allele sizes at markers 
located between markers D18S469 and D18S554, inclusive, and 

determining the most effective treatments for said genotypes of affected individuals. 
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7. A method of detecting the presence of a bipolar mood disorder 
susceptibility locus in an individual comprising: 

analyzing a sample of DNA from said individual for the presence of a 
DNA polymorphism on the long arm of chromosome 18 between D18S469 and 
D18S554, wherein said DNA polymorphism is associated with a form of bipolar 
mood disorder. 

8. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between D18S469 and D18S1161, inclusive. 

9. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between D18S1161 and D18S1 12 1 , inclusive. 

10. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between D18S1121 and D18S1009, inclusive. 

11. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between D18S1109 and D18S380, inclusive. 

12. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between and D18S380 and D18S554, inclusive. 

13. The method of claim 7, wherein said DNA polymorphism is located on the 
long arm of chromosome 18 between D18S1009 and D18S554, inclusive. 

14. The method of claim 7, wherein said analyzing further comprises: 

a. obtaining DNA samples from family members of said individual, 

b. analyzing said DNA samples from family members for the presence of 
said DNA polymorphism, and 

c. correlating the presence or absence of the DNA polymorphism with 
a phenotypic diagnosis of bipolar mood disorder for said individual and for said 
family members. 
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15. A method for detecting the presence of a DNA polymorphism linked to a 
gene associated with bipolar mood disorder in an individual comprising: 

a. typing blood relatives of said individual for a DNA polymorphism 
located within a region of chromosome 18, wherein said region is located between 
D18S469 and D18S554, inclusive, 

b. analyzing a DNA sample from said individual for the presence of 
said DNA polymorphism. 

16. A method of genetically diagnosing bipolar mood disorder in an individual 
comprising: 

a. obtaining a DNA sample from said individual, 

b. analyzing said DNA sample for the presence of a DNA 
polymorphism associated with bipolar mood disorder, wherein said DNA 
polymorphism is located within a region of chromosome 18, wherein said region is 
located between D18S1121 and D18S554, inclusive. 

17. A method of confirming a phenotypic diagnosis of bipolar mood disorder in 
an individual comprising: 

a. obtaining a DNA sample from said individual, 

b. analyzing said DNA sample for the presence of a DNA 
polymorphism associated with bipolar mood disorder, wherein said DNA 
polymorphism is located within a region of chromosome 18, wherein said region is 
located between D18S1121 and D18S554, inclusive. 

18. The method of claim 17, wherein said individual has Spanish or 
Amerindian ancestry. 
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19. A method of classifying subtypes of bipolar mood disorder comprising: 

a. identifying one or more DNA polymorphisms located within a 
region of chromosome 18, wherein said region is located between D18S1121 and 
D18S554, inclusive; and 

b. analyzing DNA samples from individuals phenotypically diagnosed 
with bipolar mood disorder for the presence or absence of one of more of said 
DNA polymorphisms. 

20. A method of treating an individual diagnosed with bipolar mood disorder 
comprising: 

a. identifying one or more DNA polymorphisms located within a 
region of chromosome 18, wherein said region is located between D18S1 121 and 
D18S554, inclusive; and 

b. analyzing DNA samples from individuals phenotypically diagnosed 
with bipolar mood disorder for the presence or absence of one of more of said 
DNA polymorphisms, and 

c. selecting a treatment plan that is most effective for individuals 
having a particular genotype within said region of chromosome 18. 

21. An isolated polynucleotide capable of selectively hybridizing with a DNA 
sample from an individual phenotypically diagnosed with severe bipolar mood 
disorder, wherein said polynucleotide does not selectively hybridize with a DNA 
sample from an individual not affected by severe bipolar mood disorder, wherein 
said isolated polynucleotide selectively hybridizes with a complementary 
polynucleotide within a region of chromosome 18, wherein said region is located 
between D18S1121 and D18S554, inclusive. 

22. The isolated polynucleotide of claim 21, wherein said complementary 
polynucleotide is within a region of chromosome 18, between D18S1 121 and 
D18S554, inclusive. 
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